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ABSTRACT 

We present new results on the evolution of the cosmic star formation rate as a function of stel¬ 
lar mass in the SXDS-UDS field. We make use of narrow-band selected emission line galaxies 
in four redshift slices between z= 1.46 and z = 0.63, and compute stellar masses by fitting a 
series of templates to recreate each galaxy’s star formation history. We determine mass-binned 
luminosity functions in each redshift slice, and derive the star formation rate density (Psfr) 
as a function of mass using the [O Ill] or [O ll] emission lines. We calculate dust extinction 
and metallicity as a function of stellar mass, and investigate the effect of these corrections on 
the shape of the overall Psfr(M). We find that both these corrections are crucial for deter¬ 
mining the shape of the Psfr(M), and its evolution with redshift. The fully corrected Psfr(M) 
is a relatively flat distribution, with the normalisation moving towards lower values of Psfr 
with increasing cosmic time/decreasing redshift, and requiring star formation to be truncated 
across all masses studied here. The peak of Psfr(M) is found in the 1O**''^<M0<1O'^ ** mass 
bin at z = 1.46. In the lower redshift slices the location of the peak is less certain, however 
low mass galaxies in the range 10^ *’<Mq<10^ *' play an important part in the overall Psfr(M) 
out to at least z 1.2. 

Key words: cosmology:observations - surveys - galaxies:evolution - galaxies:formation - 
galaxies:high-redshift - galaxiesduminosity functions. 


1 INTRODUCTION 

The evolution of cosmic star formation is well-studied, and cru¬ 
cial to our understanding of the Universe as a whole. It is now 
well-established that the peak of star formation activity lies beyond 
z ~ 1 (e.g. Hopkins & Beacom 2006) and that there has been a 
steep decline in this activity to the present day (Lilly et al. 1996, 
Le Floc’h et al. 2005, Prescott et al. 2009, Sobral et al. 2013, Drake 
et al. 2013). Determining the cosmic star formation rate density as a 
function of stellar mass however, Pspr(M), is somewhat more diffi¬ 
cult, and the specifics of this aspect of observational cosmology are 
still obscure. Reports of a shift in the masses of star-forming galax¬ 
ies across redshift are now commonplace, however the meaning of 
this is complex. Stars in massive galaxies are known to have formed 
earlier and on a shorter timescale than those formed in lower mass 
systems (Stanford et al. 1998, Heavens et al. 2004, Panter et al. 
2007), and while this switch-off of star formation in massive galax¬ 


ies is readily observed (e.g. Brodwin et al. 2013) and attributed 
to the feedback of an active galactic nucleus (AGN) truncating star 
formation, the behaviour of the lower mass star-forming population 
is less clear. For instance, the shift of the primary location of star 
formation activity from higher to lower mass galaxies with increas¬ 
ing time is clearly seen when considering specific star formation 
rates (sSFRs; star formation rates per unit stellar mass e.g. Noeske 
et al. 2007), however, this places only a weak constraint on the 
masses of objects dominating the overall Psfr at each z- 

In addition to AGN feedback in high mass galaxies, super¬ 
novae are known to expel gas from lower mass systems preventing 
the continuation of star formation (Oppenheimer & Dave 2008), 
and recent work has suggested that ~ 3 percent of star forming 
galaxies at z ~ 1 may have their star formation truncated through 
outflows of gas (e.g. Yuma et al. 2013). Determining the role of 
low-mass galaxies and establishing the true shape of the overall 
Psfr(M) is therefore the first step towards uncovering the relative 
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importance of the physical processes at work in truncating star for¬ 
mation. 

The Redshift One LDSS3 Emission Line Survey (ROLES; 
Davies et al. 2009) was established to address the role of low-mass 
galaxies in the overall PsFR a function of mass and selected a 
sample of dwarf galaxies 1O^'^<M0<1O®'^ at z ~ 1. Gilbank et al. 
(2010) used [O II] emission to estimate star formation rates (SFRs) 
for these galaxies and concluded that the shape of Psfr(M) does 
not evolve with redshift between z ~ 1 and the present day, a result 
since corroborated by Peng et al. (2010) and Sobral et al. (2014) 
who have both reported little change in the mass function of star¬ 
forming galaxies since z ~ 1. 

The strong [O II] emission line doublet (3726A and 3729A) 
used in Gilbank et al. (2010) provides a particularly useful tool to 
trace star formation at redshifts >1 where Ha is shifted out of the 
optical window, and yet constraints on SFR via [O II] emission are 
estimated with a number of caveats. The traditional conversion be¬ 
tween observed L[0 II] and SFR (Kennicutt 1998) determines the 
average [O Il]/Ha ratio for a sample of 90 local star-forming galax¬ 
ies, and uses this to convert the apparent [O II] luminosity to an 
Ha-derived estimate of the SFR. Applying this relation to different 
samples of galaxies however, particularly if these are split by stel¬ 
lar mass, renders the relationship invalid as the average reddening 
from Ha to [O II] of samples at different masses may be drastically 
different to local star-forming galaxies. 

Extinction is a frequently occurring obstacle for extragalactic 
astronomy, and many studies have attempted to quantify the ef¬ 
fects and dependence of extinction on other physical parameters, 
e.g. stellar mass, metallicity or SFR (Heckman et al. 1998, Hopkins 
et al. 2001, Brinchmann et al. 2004, Asari et al. 2007, Garn et al. 
2010a). Using the Balmer decrement, Garn & Best (2010) found 
that for a star-forming galaxy, extinction correlates most strongly 
with stellar mass, and propose a new relationship to describe this 
relationship calibrated for Ha luminosity but which can be extrap¬ 
olated to other wavelengths. 

As an additional consideration when using a metal line to 
determine SFR, one must account for the effect of the gas-phase 
metallicity on the strength of the lines. The [O III] and [O II] for¬ 
bidden lines for instance are sensitive to metal abundances in ad¬ 
dition to the temperature of the gas (Jansen et al. 2001, Chariot 
et al. 2002). In previous work Gilbank et al. (2010) for example 
applied an empirical correction to simultaneously account for the 
effects of both extinction and metallicity on the [O II] indicator. 
As an alternative however, Kewley, Geller & Jansen (2004) deter¬ 
mine a conversion between intrinsic L[0 II] and SFR which is thus 
independent of the reddening between the wavelengths of Ha and 
[O II], and also includes an optional correction for the metallicity of 
the gas. This relationship can then be applied to the intrinsic [O II] 
luminosities (i.e. those corrected for extinction at the wavelength 
of [O II]) of any sample of galaxies. 

Understanding the evolution of Psfr(M) requires a dataset 
with a number of attributes. In addition to the requirement of a 
reliable SFR indicator, the data must be deep to probe low SFR ob¬ 
jects, as well as sampling a wide range of masses. This problem is 
well suited to a narrow-band-selected sample where objects are de¬ 
tected via line emission across a wide range of masses and are close 
to SFR-limited (e.g. Sobral et al. (2014) who select star-forming 
galaxies based on their Ha emission to examine the evolution of 
the SFR and mass functions between z = 2.23 and z = 0.40). 

In this paper, we utilise the ability of narrow-band selection to 
detect large numbers of star-forming galaxies irrespective of their 
stellar mass, to allow a statistical evaluation of the objects’ prop¬ 


erties as a function of mass. We build on the work of Drake et al. 
(2013) by determining stellar masses for objects in the four redshift 
slices covered by deep Subaru data at z = 0.63, z = 0.83,z = 1.19 
and z = 1.46. We fit maximum likelihood luminosity functions to 
our sample using [O III] or [O II] emission in a series of mass bins 
to examine the resultant values of psFR- We incorporate a careful 
treatment of extinction and metallicity effects as a function of stel¬ 
lar mass, calibrating the dependence of SFR on metallicity down to 
10^ ** < M© < 10*'®. We aim to highlight the effect of these cor¬ 
rections on the shape of the Psfr(M) and consider the implications 
for the physical processes required to produce the fully corrected 
Psfr(M). 

This paper proceeds as follows: in Section 2 we describe in 
brief the data used for this analysis, in Section 3 we explain the 
methods used to determine stellar masses, produce mass-binned 
luminosity functions, calculate star formation rates and the correc¬ 
tions applied to allow a treatment of extinction and metallicity as 
a function of stellar mass. In Section 4 we present our results in 
the form of luminosity functions and values of Psfr psr mass bin 
before discussing the factors affecting our results and the scientific 
implications of the shape of the fully corrected Psfr(M) between 
redshifts 1.46 and 0.63. We then draw conclusions in Section 6. 

For consistency with the manner in which we calculate stellar 
masses all values of Psfr are calculated for a Chabrier IMF. An Hq 
= 70 kms^* Mpc^*, Q.m =0.3 and flyv =0.7 cosmology is assumed 
throughout, and all magnitudes are in the AB system. 


2 DATA AND SAMPLE SELECTION 

We use data from the Subaru/XMM-Newton Deep Survey (SXDS; 
Furusawa et al. 2008) and the UKIDSS Ultra Deep Survey (UDS; 
Foucaud et al. 2007). The samples discussed in this paper are 
drawn from the catalogue of Drake et al. (2013), and are se¬ 
lected using two narrow-band filters on the Subaru Telescope at 
8150 A (NB816) and 9183 A (NB921). The observations are com¬ 
plemented by 11 bands of photometry: CFHT n-band, Subaru B, 
V R, i and z bands, UKIRT J, H and K bands, and Spitzer IRAC 
coverage in channels 1 and 2 (see table 1 of Drake et al. 2013 for 
further details). The narrow-band imaging has been smoothed us¬ 
ing a Gaussian kernel, in order to match the point spread function 
(PSF) of the SXDS broad-band data (Ouchi et al. 2008, 2009). 

Objects are selected according to their narrow-band excess, 
requiring a 5 ct detection in the narrow band, a colour excess of 
at least 3 ct above the sigma-clipped median intrinsic colour, and 
of at least 3(T relative to scatter introduced through photometric 
uncertainty. 

The sample of line-emitters is cleaned of late-type stellar con¬ 
taminants via the BzK technique of Daddi et al. (2004), and pho¬ 
tometric redshifts are determined using the photometric redshift- 
fitting software “EAzY”, making full use of the 11-band photom¬ 
etry to assign objects to redshift slices. This approach results in 
photometric redshifts with a normalised median absoulte deviation 
of cTnmad = 0.026. 

Finally, the completeness of detection in the narrow-band 
sample is assessed per AmNB = 0.05 bin using randomly positioned 
fake objects. The stringent 5 ct detection limit means the effect is 
small, and only small numbers of objects need to be accounted for 
in the maximum likelihood analysis. 

Full details can be found in Drake et al. (2013) of data cover¬ 
age, sample selection, photometric redshift determination and as¬ 
signment to redshift slices (via stacked probability density distri- 
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Table 1. Objects per redshift slice 


Filter 

Redshift slice 

Emission Line 

Objects 

NB816 

0.35 < 0.63 < 0.80 

[0 III] 

999 

NB921 

0.50 < 0.83 <1.10 

[0 III] 

894 

NB816 

0.80 < 1.19 < 1.50 

[0 II] 

956 

NB921 

1.10 < 1.46 < 1.90 

[0 II] 

2203 


butions) plus further information on the maximum likelihood ap¬ 
proach to luminosity functions. Table 1 gives numbers of objects in 
each redshift slice used here. 

3 METHOD 

3.1 Stellar Masses 

We determine stellar masses following the method of Simpson et al. 
(2012). Each object’s redshift is set to the median redshift of the 
redshift slice, and the SED is fit with a series of synthetic star for¬ 
mation histories (SFHs) using the redshift-fitting sofware EAzY 
to determine the most probable scenario for the galaxy’s assem¬ 
bly. Simpson et al. (2012) design a set of 40 SFH templates based 
on instantaneous starbursts using a Chabrier IMF and the spectral 
models of Chariot and Bruzual (2007). The SFH templates range 
between very young populations, ~ 20 x 10 ^ years old, up to popu¬ 
lations 13x10® years old. Finer time resolution is used between the 
younger templates, and a few young reddened templates are also 
included. The advantage of this approach is the ability to recreate 
any SFH through combination of instantaneous starburst templates. 
As this calculation is based on 2 arcsecond aperture fluxes from the 
catalogue described in Drake et al. (2013) however, each mass must 
be scaled by the ratio of total K-band flux to aperture K-band flux. 
Here we use total K flux values from the catalogue of Grutzbauch 
et al. ( 2011 ) where available, or a 2.2 arcsecond aperture flux to 
simulate total K-band flux otherwise. 

Figure 1 shows a comparison of masses derived in this work 
and the stellar mass estimates of the same objects detected in Hart¬ 
ley et al. (2013). Not all of the star-forming galaxies detected in 
our survey are found in the K-band selected sample of Hartley et al. 
(2013), however where the samples overlap our stellar masses agree 
well with the published work. Our use of the Chariot and Bruzual 
(2007) models which incorporate a greater contribution from ther¬ 
mally pulsing (TP-) AGB stars, results in lower masses by 0.1 — 0.2 
dex than those of the Bruzual & Chariot 2003 models used in Hart¬ 
ley et al. (2013). 

3.2 Mass-Binned Maximum-Likelihood Luminosity 
Functions 

Drake et al. (2013) described our method of determining maximum 
likelihood luminosity functions for narrow-band selected samples, 
and demonstrated that this approach possesses a number of advan¬ 
tages over existing techniques. We now apply this method to mass- 
binned data at z = 0.63, z = 0.83, z = 1.19 and z = 1.46 to deter¬ 
mine luminosity functions as a function of stellar mass. A small 
but important difference between this analysis and that of Drake 
et al. (2013), is the manner in which we model the underlying star¬ 
forming galaxy population. For the original analysis the population 
of broad-band selected galaxies in a redshift slice was well fit by 
a Schechter function, however splitting the sample by stellar mass. 






Figure 1. Comparison of stellar mass estimates from this work and results 
in the literature from Hartley et al. 2013. The four panels represent the 4 
redshift slices we have used for this analysis progressing from top left to 
bottom right in increasing redshift order z = 0.63,0.83,1.19,1.46. 


this approach is obviously no longer appropriate. We find the best 
fit when this population is modelled as a Gaussian (see Figure 2), 
and so we use this to represent the broad-band magnitudes of galax¬ 
ies in each mass range from the K-band selected catalogue of Simp¬ 
son et al. (2012). For each redshift slice we fit to the observed broad 
band corresponding to the rest-frame wavelength of the NB selec¬ 
tion filter, i.e. for NB816-selected objects this is an interpolation 
of the i and z bands, and for NB921-selected objects z-band only. 
The Simpson et al. (2012) catalogue is used here in preference to 
Grutzbauch et al. (2011) or Hartley et al. (2013) for consistency in 
stellar mass estimates only. 

Each redshift slice is split into 5 — 7 mass bins (depending on 
the number of detections and range of stellar masses) of width 1 
dex in the lowest mass bin ( 1 O®'®<M 0 < 1 O* ®) and 0.5 dex for the 
remainder of the sample. 

3.3 Star Formation Rates 

To investigate the effects of different assumptions about dust and 
metallicity on the resultant SFRs, we compute pspR using a series 
of different approaches. 

To determine SFR from Lpm] we initially take the same ap¬ 
proach as in Drake et al. (2013) and apply the standard line ratios: 
Ha/Hj6=2.78 and [O IIl]/Hj3=3 (Osterbrock & Ferland 2006) to 
convert the Kennicutt (1998) relation between Lna and SFR to one 
for Ljoni] (Equation 1): 

SFR(M 0 yr-*) = 7.35 x 10-35L[o,u]E[om] (1) 

where Ljqjjjj is the observed [O III] luminosity, and Epm] repre¬ 
sents extinction at the wavelength of [O III]. The Kennicutt (1998) 
approximation of SFR for Ljon] is computed according to Equation 
2 : 

SFR(M 0 yr-') = 1.39 x 10-34 L[o,i]Eh„ (2) 
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z band 


Figure 2. Example of Gaussian fit to the distribution of broad-band mag¬ 
nitudes at z = 1.46 in the 10.0 < logiM/Mp,) < 10.5 mass bin. For these 
NB921-selected objects this is the z band magnitude, taken from the K- 
band selected catalogue of Simpson et al. (2012) 


where Ljqjjj is the observed [O II] luminosity, and Ena represents 
extinction at the wavelength of Ha. 



Figure 3. Metallicity per mass bin of width 0.5 dex. Metallicity is particu¬ 
larly low in the lower half of mass bins. Solar abundance in these units is 
8.9 marked by the dark orange dotted line. The dashed horizontal line at 8.4 
represents the limit of the metallicities modelled in the Kewley et al. (2004) 
calibration for SFR. 


extinction at the wavelength of [O II]. We apply this conversion to 
determine pspR for galaxies with a dust extinction correction as a 
function of stellar mass, using Equation 3 in conjunction with the 
prescription of Cardelli et al. (1989) to determine Ejou]- 


3.3.1 Mass- Dependent Extinction Correction 

Extinction is known to correlate with a number of physical prop¬ 
erties of galaxies (Brinchmann et al. 2004, Gam et al. 2010b) the 
dominant of these being stellar mass (Garn & Best 2010). To ad¬ 
dress the mass-dependence of dust extinction, we implement Equa¬ 
tion 3 of Garn & Best (2010) to describe Ha extinction in magni¬ 
tudes as a function of stellar mass: 

Ana =0.91+0.77x4-0.1 1 x:2-0.09x:^ (3) 

where x = logio(M/lO*®M 0 ). Using Equation 3 we determine val¬ 
ues for Ahcs per mass bin and apply the Cardelli et al. (1989) red¬ 
dening law to the set of values found for this prescription to deter¬ 
mine the corresponding values of Ajqjjjj and Ajoqj . Each fit is set 
to a constant below IO^'^Mq where the estimate is unreliable. 

For Ljoni] estimates of SFR we use Equation 1 in conjunc¬ 
tion with an assessment of the dust extinction per mass bin from 
Equation 3. For L[on] however, the estimate of SFR from Equation 
2 is affected by uncertainties due to the wavelength dependence of 
reddening between the Ha and [Oil] lines. While this is the stan¬ 
dard approach for large statistical samples, it is unsuitable for this 
dataset, binned by stellar mass and spanning a wide range of red- 
shifts. Kewley et al. (2004) re-calibrate the [Oil] indicator to be 
reddening independent, allowing for the application of SFR([OII]) 
to a wide range of galaxy samples. This is presented in Equation 4: 


3.3.2 Metallicity Correction 

Although the [O III] and [O II] forbidden lines are sensitive to the 
the metallicity of a galaxy. Figure 11 of Kewley et al. (2004) shows 
that the ratio of [O IIl]/Hj3 line luminosities is constant down to 
very low metallicities for a range of ionisation states of the gas, 
and hence we apply no metallicity correction for [O IIl]-derived 
estimates of SFR. 

[O II] emission however varies considerably with gas- 
phase metallicity (here meaning oxygen abundance defined as 
log(0/H) -I- 12), and since the traditional conversion of Ljo^j to 
SFR relies on the relationship between observed [O II] to Ha, the 
variation with metallicity introduces an error in this approach. Kew¬ 
ley et al. (2004) build on their earlier work to incorporate the effect 
of metallicity on the [O Il]/Ha ratio, and hence derive a theoretical 
calibration for SFR(L|on]) a function of L|on] and metallicity. * : 


SFR([OII],Z)t(M 0 yr-') 


7.9x 10-35 l([OII]) 
[Oil]/Ha 


(5) 


They use theoretical models to examine the dependence of the 
[O Il]/Ha ratio on metallicity for a number of different ionisation 
states of the H II gas. The relationship for intrinsic [O Il]/Ha line 
ratio on metallicity is then given by: 


[Oil ]/Ua = 0 3 - bZ + cZ^ + dZ^ 


( 6 ) 


SFR([OII])(M 0 yr-i) = (6.58+1.65) x lO-^^LfonjEpn] (4) 
where L|on] is observed [O II] luminosity, and E|oii] is the flux 


where Z= log(0/H) -(-12, and and the coefficients a, b, c, and d orig¬ 
inate from the best fit curve for metallicity abundance vs [O Il]/Ha 
for the appropriate ionisation parameter. 
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For consistency with our assumption of [O IIl]/Fi/3=3 in Drake 
et al. (2013), we adopt an ionisation parameter- q = Ax 10^. This 
is equivalent to log U = —2.875, where (/, the commonly adopted 
ionisation parameter, is defined as U = q/c. Our choice is en¬ 
tirely consistent with recent detailed studies e.g. Nakajima & Ouchi 
(2014) who found ^ = 1 x 10^ — lx 10*, and also falls within the 
“normal” range of ionisation parameters found in Kewley et al. 
(2004), and Dopita et al. (2000): ^ = 1 x 10^ — 8 x 10^ cm s^'. 
For ^ = 4 X 10^ a = -1432.67, b = 470.545, c = -51.2139 and 
d= 1.84750. 

To derive metallicities for our sample we employ the empiri¬ 
cal model of Savaglio et al. (2005) for the evolution of the mass- 
metallicity relation to high redshift. They derive an expression for 
the metallicity of a galaxy of a given stellar mass, at t/f, the age of 
the Universe at that redshift for the assumed cosmology, as: 

12-Flog(0/H) = -7.5903 + 2.5315 logM 
-0.09649 log2M + 5.1733 log fn -0.3944 log^ 

—0.4030 log tH log M. 

(7) 

We use Equation 7 to determine an average metallicity per 
mass bin, for each redshift slice. The resultant values can be seen in 
Figure 3, and reach substantially sub-solar values. The metallicity 
of the Sun is shown on this plot as the dotted orange line at Z=8.9. 

For mass bins of very low metallicity (< 8.2, about half of 
our sample) the model fit to [O Il]/Ha from Kewley et al. (2004) 
breaks down, and extrapolating the fit results in non-physical val¬ 
ues of [O Il]/Ha. The limit of the Kewley et al. (2004) model is 
shown as the dashed black line in Figure 3. Kewley et al. (2004) 
use Pegase stellar population models (Fioc & Rocca-Volmerange 
1997) with the Starburst99 code (Leitherer et al. 1999) to deter¬ 
mine an ionising spectrum to simulate values of [O Il]/Ha, using 
the Mappings hi radiative transfer code (e.g. Sutherland & Dopita 
1993) for different values of metallicity and ionisation parameter, 
q. 

To determine values of [O Il]/Fla we therefore follow the 
method of Kewley et al. (2004) but use the CLOUDY (Ferland 
et al. 2013) radiative transfer code to determine [O Il]/Ha at each 
metallicity. We assume a total H density of 10^'^ cm^* which cor¬ 
responds to the electron density 350 cm^* used in Kewley et al. 
(2004). The resultant values of [O Il]/Ha agree to ~ 10% over the 
corresponding metallicity range. 

The [O Il]/Ha ratios are well fit in log space by Equation 8 : 

log([01I]/Ho:) =fl + foZ + cZ^ + dZ* (8) 

with best fitting coefficients a = 76.94, b = —33.86, c = 4.77, and 
d = -0.22. 

We use Equation 5 in conjunction with our own estimates of 
[O Il]/Ha to determine fully corrected values of pspR based on 
[O II] luminosity. 


^ where q is the maximum velocity of the ionisation front driven by the 
local radiation field (Kewley et al. 2004) 


4 RESULTS 

4.1 Luminosity Functions 

The mass-binned maximum-likelihood luminosity functions for 
objects at the weighted-mean redshift of emitters in a particular red¬ 
shift slice are presented in Figure 4 and Table AL Volumes quoted 
in the table are representative of the volume where narrow-band fil¬ 
ter transmission is greater than 0.5. Values of L* are corrected for 
Galactic extinction and aperture effects. The data in Figure 4 are 
presented as in Drake et al. (2013), where the solid coloured line 
gives the maximum likelihood fit, and the data points on the plot 
are binned arbitrarily for presentation. 

For each luminosity function the maximum likelihood analy¬ 
sis determines the most likely values of (j)*, L* and a, allowing all 
three parameters to vary. In order to produce a well-constrained fit 
however, the data must probe significantly below L* in the particu¬ 
lar redshift/mass-bin combination. For some luminosity functions, 
the depth of the data coupled with a small number of emitters in 
that mass bin, result in poorly-constrained values of a. Where this 
is the case, we follow the approach taken in Drake et al. (2013), 
and constrain the luminosity function with a Gaussian prior on the 
faint end slope, set to the value of a for the closest mass bin in 
that redshift slice. A point of note is that in each case the sum of 
the mass-binned luminosity functions (denoted by the dashed line) 
including those with a Gaussian prior, is in good agreement with 
the best-fit luminosity function for the entire redshift slice (given 
by the solid black line). The small deviations between the two lines 
seen at the bright or faint end of the LFs can be explained by ob¬ 
jects in the wings of the mass distributions where there were too 
few sources to fit a LF. 

4.2 psFR as a Function of Stellar Mass 

Figure 5 presents values of Psfr(M) evaluated considering an in¬ 
tegration of the data to two different limits in the left and right 
hand panels, and incorporating various different approaches to the 
treatment of dust extinction and metallicity. The method used for 
each row of panels is detailed below and the corresponding values 
are summarised in Table Bl. For consistency with the manner in 
which we calculate our stellar masses, the PsFR values quoted and 
plotted have been converted to a Chabrier IMF. 

Initially, psFR is computed using the Kennicutt (1998) relation 
and applying 1 magnitude of extinction at Ha (upper panels Fig¬ 
ure 5). In the central panels of Figure 5, we determine extinction 
as a function of stellar mass according to the prescription of Gam 
& Best (2010) using Cardelli et al. (1989) to infer values at [O III] 
and [O II], and use Equation 4 of Kewley et al. (2004) to derive 
SFR([0 II]). Einally, in the lower two panels we use our estimates 
of log(0/H) -F 12 to incorporate a metallicity correction in conjunc¬ 
tion with Equation 5 of Kewley et al. (2004) to determine a fully 
corrected Psfr(M). 

The appropriate limits of integration to determine PsFR psr 
mass bin are complicated for this analysis, since the sensitivity in 
M 0 yr^* varies considerably due to the use of two different emis¬ 
sion line indicators and the range of redshift studied. Here we 
evaluate each luminosity function in two different ways; to the 
limit of the data in each redshift slice (left-hand panels Figure 5) 
and to the same limiting SFR in each redshift slice (right-hand 
panels Figure 5). In the former approach, the limiting SFRs are 
0.08,0.12,0.55 and L42M0yr-' at z = 0.63,0.83,1.19 and 1.46 
respectively, producing sensible results for most mass bins. Inter¬ 
grating to the same limiting SFR for comparison across redshift 
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Figure 4. Best-fitting luminosity functions for redshift slices: z = 0.63 between 7.25 and 9.25 dex, z = 0.83 ( 7.75 to 10.25 dex), z = 1-19 (7.75 to 10.75 dex), 
z = 1.46 (7.75 to 11.25 dex). The upper SFR scale on these plots gives an indication of SFR values assuming a constant scaling from luminosity for these LFs. 
Neither these numbers nor the luminosities have been corrected for dust extinction. 


slices however presents more of a challenge. Selections made in 
the z = 0.63 and z = 0.83 redshift slices can be integrated to low 
SFRs ~ 0. lM©yr^ *, however limits which produce sensible results 
for these redshift slices require a large extrapolation at z = 1.19 
and z = 1.46 and consequently introduce a far greater uncertainty 
on luminosity functions that probe little below L*. Likewise, limit¬ 
ing the integration to values where the luminosity function is well- 
constrained in these redshift slices results in the loss of a large por¬ 
tion of the Psfr(M) in the lower redshift slices (additionally there 
is a dramatic drop in the 1 O®'^<M 0 < IqIO.O at z = 0.63 where 
the whole luminosity function is poorly constrained) and ultimately 
the best choice for a constant limiting SFR is IMQyr^ *. We include 
this analysis to provide a consistent evaluation across redshift, how¬ 
ever we prefer values of Psfr(M) based on the limit of the data as 
being representative of the total SFR. 


4.2.1 Errors on pspj^(M) 


Errors on PsFR ars first computed according to the maximum and 
minimum values that arise from the IcT deviation of (/>*, L* and a 
when all three parameters are allowed to vary in the maximum- 
likelihood analysis. Additionally we quantify the uncertainty in 
PSFR introduced through the small number of objects making up 
some mass bins by incorporating an additional Poissonian error. 
The initial 1 a error is combined with the fractional error on psFR 
according to Gehrels (1986) given the number of objects making up 
the luminosity function. Errors quoted in Table B1 and shown on 
Eigure 5 represent the combination of these two errors in quadra¬ 
ture. 
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Figure 5. psFR as a function of mass for redshift slices at z = 0.63, z — 0.83, z— 1-19 and z = 1 -46. psFR is estimated via Kennicutt (1998) and 1 mag extinction 
at Ha (top two panels), Kewley et al. (2004) and dust extinction as a function of stellar mass (central two panels), and the fully connected psFR including dust 
an metallicity corrections (lower two panels). The left-hand column of figures present results for an integration to the limit of the data in each redshift slice, 
and the right hand column of figures presents an integration to lM 0 yr“^ 
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5 DISCUSSION 

5.1 The Effect of Dust and Metallicity Corrections on 

Psfr(M) 

Applying 1 magnitude of extinction at Ha (top two panels of Figure 
5) shows the peak of Psfr(M) lying in the 1 O®'®<M 0 < 1 O®-^ bin at 
z = 1.46 andz = 1.19, shifting to 10 *-^<Mq< 10®'® at z = 0.83, and 
z = 0.63, a feature which is preserved for both limits of integration. 
The most notable difference between the two approaches is the far 
lower values of Psfr(M) in the two lower redshift slices due to the 
lM 0 yr^' limit probing little below L*. 

Incorporating dust corrections as a function of stellar mass 
(central two panels) produces a dramatic change in the shape of the 
Psfr(M) for the redshift z = 1.46 slice, raising the high mass end 
(>1O*®'*’M0) significantly and lowering values of Psfr(M) for all 
masses < 1 O*®'®M 0 . Similarly the high-mass portions of Psfr(M) 
in the z = 1.19 and z = 0.83 redshift slices are lifted, however in 
the z = 0.63 slice, little change is seen due to the lack of massive 
objects. 

The lower two panels of Figure 5 present our best es¬ 
timates of Psfr(M) incorporating dust and metallicity correc¬ 
tions. Little change is seen across the central region of the func¬ 
tion (10**''<M 0 < 1 O''*'') however a small rise is seen in masses 
> 1 O*'’''’M 0 , and a large increase below 1O* *'M 0 . This low mass in¬ 
crease only affects the z = 1.19 slice since galaxies below 1 O^'*’M 0 
do not enter the z = 1.46 selection, and the low mass objects at 
z = 0.63 and z = 0.83 have their SFRs computed via [O III] lumi¬ 
nosity which requires no metallicity correction. 

Figure 6 presents in a simple format the way each correction 
applied to the Psfr(M) affects the shape of the function. In each 
panel the black shaded area and dashed line highlight the shape of 
the Psfr(M) using the traditional assumptions about [O Il]/Ha and 
the converison of Kennicutt (1998). The dashed orange line then 
gives the Psfr(M) corrected only for dust extinction as a function 
of stellar mass, and the red dashed line gives the final Psfr(M) 
corrected for extinction and metallicity (where appropriate) as a 
function of stellar mass for each redshift slice. 


5.2 Evolution of Psfr(M) 

We favour the lower left-hand panel of Figure 5 as representative 
of the fully-corrected Psfr(M). The peak contribution to Psfr at 
z= 1.46 comes from galaxies in the mass range 1 O"''^<M 0 < 1 O" '', 
however in the lower redshift slices it is less clear where the peak of 
the function lies. The objects detected in the lower redshift slices do 
not reach masses as high as those at z = 1.46 and so it is impossible 
to say if the psFR is still rising at these masses. Figure 6 highlights 
how dust extinction as a function of mass kicks in at the highest 
masses which would act to raise the Psfr considerably at the high 
masses which are missing from the lower two redshift slices. 

The shape of the fully corrected Psfr(M) is remarkably flat 
in comparison to the raw function, and shifts to lower normalisa¬ 
tion with increasing time/decreasing redshift. A similar effect is 
seen in both Gilbank et al. (2010) and Sobral et al. (2014) who 
conclude that the shape of the Psfr ^s a function of mass shows 
very little evolution with redshift. One thing we can note with this 
particular selection of galaxies is that galaxies of very low mass 
(1O^'''<M0<1O^ '') play an important part in the contribution to the 
overall Psfr out to at least z = 1.19. 



Figure 6. The shape of the pspR as estimated via Kennicutt (1998) and 1 
mag extinction at Ha (black shaded region plus black dashed line), Kewley 
et al. (2004) and dust extinction as a function of stellar mass (orange shaded 
region and orange dashed line), and the fully con'ected Psfr including dust 
an metallicity corrections (red shaded region and red dashed line). 


5.3 The Contribution of Different Masses to Psfr 

Considering galaxies in four mass bins where we have 
data points for all four redshift slices: 1 O^'’<M 0 < 1 O*-^ 
(cyan), 1 O*'^<M 0 < 1 O®'*’ (green), 1 O®''<M 0 < 1 O®'^ (yellow), 
1 O^'^<M 0 < 1 O*'''’ (orange). Figure 7 examines the contributions 
from galaxies of different masses to the overall Psfr- The dotted 
line and green Z symbols represent the sum of these four bins, 
and coloured points represent the contributions to Psfr from each 
of these mass bins. Broadly speaking the relative contributions 
from the different mass bins decline simultaneously with the 
value of Psfr itself. A notable feature however is the stronger 
decline in contribution from galaxies in the 1 O®'^'’<M 0 < 1 O**'''’ bin 
relative to the lower mass bins. This is in contrast to the result 
of Sobral et al. (2014) where Psfr of galaxies in all three mass 
bins presented in their figure 8 decline contemporaneously with 
the overall Psfr- There are likely a number of reasons which lead 
to this discrepancy. Sobral et al. (2014) benefits from a consis¬ 
tently selected sample of Ha emitters across redshifts z = 0.40 
- z = 2.23 meaning SFR is calculated homogeneously across 
the entire redshift range. In contrast we have estimated Psfr(M) 
using [O III] (z = 0.63,z = 0.83) or [O ii] (z = 1.19,z = 1.46) 
emission, indicators which are known to have a higher and less 
well understood contamination from AGN (see Section 5.4). 

Additionally, slightly different (larger) mass bins are used for 
the Sobral et al. (2014) analysis, plus the larger range in redshift 
studied may allow the small fluctuations seen across the four red¬ 
shift slices studied here to be smoothed out. 
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Figure 7. Contribution of galaxies of different masses to the overall psFR 
in the four redshift slices studied here. Dot-dashed lines and points there¬ 
on represent the overall psFR taken from Drake et al (2013) converted to 
a Chabrier IMF. The dashed line and rho symbols present the sum of all 
mass bins in this study, and the dotted line and capital Sigma symbols 
show the sum of the four mass bins examined here. 1O* ®<M0<1O*'^ (cyan), 
1O*'^<M0<1O‘’“ (green), 1O® °<M0<1O‘’'5 (yellow), 10 ‘’'^<Mq< 10‘““ (or¬ 
ange) 


5.4 The Effect of AGN Contamination on Psfr 

We can estimate the level of AGN-contamination by quantifying X- 
ray sources in the field that are coincident with the positions of our 
emission-line-selected sample. Using the catalogue of Ueda et al. 
(2008) and matching within a 3 arcsecond radius, we find small 
numbers of interlopers in each redshift slice: 7, 1, 0 and 9 objects 
are found in redshift slices z = 0.63, z = 0.83, z = 1.19 andz = 1.46 
respectively, equating to 1 percent (to the nearest percentage point) 
in the lowest redshift slice, and less than this in all other bins. 

We take the same approach as in Drake et al. (2013) and 
choose not to actively correct for AGN contamination, since those 
objects with AGN activity will undoubtedly be associated with 
some star formation, thus leading to a possible over correction. 

We should however consider the consequences of possible 
AGN contamination. If, as many authors have suggested, AGN con¬ 
tamination rises towards higher redshift (e.g. 10 percent at z < 1, 15 
percent at z > 1 Garn et al. 2010a, Sobral et al. 2013) the effect on 
Psfr would be to lower estimates above z = 1 by a larger amount 
than estimates at z < 1, effectively moving the four coloured con¬ 
tours closer together, meaning there has been less evolution of the 
Psfr than implied in Figure 5. 

A second consideration is the likelihood of the higher mass 
bins being more highly contaminated by AGN than the lower mass 
bins. To place some constraint on whether this may be the case, we 
chose to compare the ratio of line flux with flux at restframe 2500 
A. While the measurements show a correlation, the highest mass 
objects show the largest discrepancy, with line fluxes that would 
produce a larger SFR than that from the restframe 2500 A flux. 
This gives an indication that the highest mass bins could be over¬ 
estimating Psfr effectively boosting the peak of Psfr at the high 
mass end. If this is the case then the shape of the Psfr(M) is in¬ 
fact flatter than depicted in Fig 5. The physical implications of this 


scenario are that the decline of star formation activity is even less 
driven by the switch off of star formation in massive galaxies, and 
governed by a process (or processes) acting across the entire mass 
range. 


6 CONCLUSIONS 

We have applied the method described in Drake et al. (2013) to 
determine mass-binned luminosity functions in four redshift slices 
between z = 1.46 and z = 0.63. We have investigated the effect of 
extinction and metallicity corrections as a function of stellar mass, 
and examined the role of galaxies of different masses in their con¬ 
tribution to the overall Psfr- Our main conclusions can be broadly 
summarised as follows: 

• The correct use of dust and metallicity corrections as a function of 
stellar mass is essential to establishing the true shape of Psfr(M) 
and its evolution with redshift. 

• The normalisation of the fully corrected Psfr(M) shifts to lower 
values with increasing cosmic time (decreasing redshift). 

• The peak in Psfr(M) is seen in the 1O*°'^<M0<1O*^'*’ mass bin 
at z = 1.46. In the lower redshift slices the location of the peak is 
less certain. 

• Low mass galaxies reaching 10^ *'<Mq< 10^'® play an important 
part in the Psfr(M). 

• The contribution to the overall Psfr from galaxies across the 
1O^°<M0<1O®'5 mass range is roughly constant between z = 1.46 
and z = 0.63. The contribution from the 1O^-^<M0<1O*°° bin 
drops more noticeably. 

In conclusion our results paint a picture in which the decline of 
cosmic star formation activity since z ~ 1 is the result of the decline 
in SFR across a broad range of masses, i.e. it can not be attributed 
simply to the switch-off of the most massive galaxies. This brings 
to the forefront the question of the physical processes required to 
execute such a decline, and the relative importance of the different 
quenching mechanisms at work in truncating star formation. 
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Table Al. Luminosity functions in mass bins for four redshift slices. Values presented here are raw results from maximum likelihood fitting. These numbers 
have not been corrected for aperture effects or dust extinction. 
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Table Bl. Values of pspR in mass bins for the four redshift slices studied here. Values of pspR nre presented first with 1 magnitude of extinction at Ha, then 
incorporating a treatment of extinction as a function of stellar mass, and finally our best estimates of pspR incorporating both dust and metallicity corrections 
(where appropriate) as a function of stellar mass. 
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:yMin p2>_2 04E-05 

^^_2.53E-05 

j.^tiP ^J_2.04E-05 

8.25 

i.^un uo_9Qg£_Q4 

4.78E-04l’;‘i:“ 

y.OUt, ‘J^_i.i9E-04 

Z.4ZC, ^^_i,i9E-04 

y.OUt ^^_i,i9E-04 

Z.4Zt, U^_i.i9E-04 

8.75 

Q f^'3+2.21E—04 

•^•^■^^"^•^-1.59E-03 

2.36E-03ll;^i:“ 

1 Ai8P a'i+ 8.63E—05 
i.08b-U3_5 ■7Q£_Q4 

i.i^P ^^_5.70E-04 

1 ^SP a'1+8.63E—05 

1.68b-03_g 2Qg_Q4 

i.i^P ^■2_5 20£_Q4 

9.25 

2.00E-03;|i=:«^ 

1 S7P A'I + I'99E—04 
1.5/b-U3_g 94£_Q4 

^■7_6.54E-04 

1 Q^p a'2 + 1'45E—04 
1.3/b-U3_g53£_Q4 

1 46E-03^*’^^^“®‘* 

l.^DP ^J_6,54E-04 

1 ':j7“C' AQ + 1.45E—04 
i.3/b-U3_gg3£_Q4 

9.75 

1.68E-03+2;‘«:“ 

1 2nF-n3+^'^^^“^ 

i.zup ^•^-5,92E-04 

1 98E-03“'“^'^^^“^‘^ 

l.yoc, ^■^-7.57E-04 

i.Hip ^^_6.98E-04 

i.yan ^j_7,57e-04 

1.7iP ^7_g9g£_Q4 

10.25 

6.68E-04+”|:04 

Z.lUC- ^^_i,04E-04 

1 31F-03+^-^^^“^^ 

i.jiP ^■J_2.91E-04 

^.iip ^‘t_2 Q4E_04 

i.jip ^^_2,91E-04 

H-iiP p‘+_2.{)4E-04 



Redshift 

z = 1.19 VoIume= 19.06 x 10“* Mpc^ 



Log mass 

1 mag Ha & Kennicutt 98 

Extinction(mass) & Kewley et al. 04 

Fully corrected: Extinction(mass) + Z 

(bin centre) 

PsFRioni (M© yr^‘ Mpc-3) 

PsFR|o,i| (M© yr^‘ Mpc-3) 

PsFR|o„| (M© yr-‘ Mpc-3) 


(>0.55) 

(>1.00) 

(>0.55) 

(>1.00) 

(>0.55) 

(>1.00) 

7.50 

H.ojn vJ‘t_2.28E-04 

O.UOP ^‘^-l,47E-04 

1.90E-04^|«i:“ 

1.24E-04;|2®:05 

1 71P A'1+6.52E—05 

1.21b-U3_j j2E_04 

i.yin '-'‘t_7.23E-05 

8.25 

1 AT+2.18E—04 

2.U/b-Ui_9 g-7£_Q4 

J1E_04 

8.45E-04^|^«i:“ 

6.08E-04i|«i:“ 

1 77P A'1+I-07E—04 
L.i /e.-UJ_ 4 g4£_04 

1 77P aq+^.SOE— 05 
1.2/b-U3_3 49£_Q4 

8.75 

/.‘+ZIZ, PJ_3.55E_03 

O.jyn ^•J-3,15E-03 

3.03E-03^|;«i:«^ 

2.69E-03;i:g|:<>4 

o.yjn 74£_03 

Q Cip rv '2 + 1.82E— 04 
3.21b-U3_j g4g_Q3 

9.25 

1 /-\o+5.35E—04 

/./jn-uo_3 ■7Q£_Q3 

7 Z7P ^t+ 5.20E—04 
7.67b-U3_3 ggg_Q3 

3.17E-03^?;‘®:U^ 

T 1 Tc a'i+2.12E— 04 
3.13b-U3_j 5Q£_Q3 

j.zyp ^J-1,81E-03 

3.26E-03^2;35E-04 

9.75 

PJ_Eg4E-03 

T Top AT+4.63E—04 

3. /oE-U3_j gjg_Q3 

2.61E-03^«5|:«4 

0 ^7P A'a+3.14E—04 
ln-VJ_^ 22,E-Q3 

7 zzp a'i+3.90E—04 

2.66b-03_j g0£_03 

2.62E-03^”^i:«^ 

10.25 

L.:)Ln ^J_^21E-04 

i.z/p ^■J_6 ,ioE-04 

i./zp ^■J_g.25E-04 

1.68E-03«;2|:« 

i.7ZP uj_99jg_04 

1 87P 04 

1.8/b-U3_9^g_Q4 

10.75 

V.ilP VJH_4 42 e_04 

O.yzn ^^_4,3 ie-04 

z.^+iP ^■J_i.i7E-03 

2.37E-03;?;|“:«^ 

3 05F-03+*-^^^“®^ 

7.U0P ^^-i,41E-03 

2 99F-03+^-^^^“*^'^ 

^.yyn 32g_Q3 



Redshift 

z = 1.46 Volume= 23.06 x 10“* Mpc^ 



Log mass 

1 mag Ha & Kennicutt 98 

Extinction(mass) & Kewley et al. 04 

Fully corrected: Extinction(mass) + Z 

(bin centre) 

PsFRioni (M© yr^' Mpc-**) 

PsFRioni (M© yr^‘ Mpc^3) 

PsFR|o„| (M© yr-‘ Mpc-3) 


(>1.42) 

(>1.00) 

(>1.42) 

(>1.00) 

(>1.42) 

(>1.00) 

8.25 

o Zip ^o+3.18E—04 

3.0lb-U3_| 72E-03 

5.37E-03li;^i:“^ 

'-'■J_7Q3£_Q4 

z.iyn PJ_i,o5E-03 

3.68E-03l^;“i:0^ 

0.46t, 26E-03 

8.75 

1 IIP AO+5.54E—04 

l-llb-O^-5.3OE-03 

1.30E-02;’|iE:M 

‘t.jjp ^■^-2.16E-03 

c Qop; AO+3.23E—04 
2.32b-U3_2 55E_o3 

6.57E-03l|^it“ 

7.71E-03^3;88i:W 

9.25 

2.25E-02+i;gi:“ 

2.45E-02ll;^i:“ 

y.L in, ^■^-4.45E-03 

1 HAP 

1.00b-02_3 Q3g_Q3 

y-y-^n 34g_03 

I.08E-02+’;«i:«^ 

9.75 

1.86E-02+‘;‘ii:“ 

1.88E-02l>;‘i:“ 

1 7^P A'7+7.54E—04 
1.2ob-U2_gQ3£_Q3 

1.27E-02+^:™:0^ 

1 7SP A9+9.06E—04 

I.z8b-Uz_2 23E-03 

i.z^p ^^_7.30E-03 

10.25 

1.33E-02+|^®:“ 

i.H-iP ^^_6,73E-03 

i./jp ^^_g.34E-03 

1.85E-02+^i3|:03 

i.oyn 00E-O2 

Z.UiP ^^_i.06E-02 

10.75 

1.28E-02+^;0>i:“ 

1 34F-n2+^‘^^“®^ 

L.j^E. ^^-6,39E-03 

3.40E-02^2;|i:“ 

3.54E-02+2|®:<>i 

H.i in ^^_i,95E-02 

H.JHP p^_2.03E-02 

11.25 

3.64E-03+“>i:“ 

^^-l,74E-03 

1.78E-02«.i|:03 

1.78E-02«;21E-03 

z.^-^p uz_j 02E_02 

z:.77P ^4i_j_Q2£_02 
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